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ABSTRACT 

Recently, Fermi-LAT detected GeV emission during the X-ray flaring activity in GRB 100728 A. We 
study various scenarios for its origin. The hard spectrum of the GeV emission favors the external inverse- 
Compton origin in which X-ray flare photons are up-scattered by relativistic electrons in the external 
forward shock. This external inverse-Compton scenario, with anisotropic scattering effect taken into 
account, can reproduce the temporal and spectral properties of the GeV emission in GRB 100728 A. 

Subject headings: gamma-ray bursts: individual (GRB 100728 A) — radiation mechanisms: non-thermal 



1. Introduction 

One of the key discoveries by the Swift satellite is 
the presence of X-ray flares during the early afterglow 
phase in a large fraction of gamma-ray bursts (GRBs) 
(see, e.g.. Burrows et al. 2005; Zhang et al. 2006; 
Nousek et al. 2006). The rapid rise and decay behav- 
ior of X-ray flares suggests that they are caused by in- 
ternal dissipation of energy due to late central-engine 
activity. It was predicted that when the inner flare 
photons pass through the forward shocks, they will 
be up-scattered by forward shock electrons, produc- 
ing a GeV flare (Wang et al. 2006; Gafli & Piro 2007; 
Fan et al. 2008). The angular dispersion effect at the 
forward shock front will wash out any shorter tempo- 
ral structure (Beloborodov 2005), so the high-energy 
inverse-Compton (IC) emission has a much smoother 
temporal structure determined by the forward shock 
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dynamical time. The incoming X-ray flare photons 
are anisotropic as seen by the isotropically distributed 
electrons in the forward shock, so the scatterings be- 
tween the flare photons and electrons in the forward 
shock are anisotropic. This effect may decrease the 
IC emission in the 1 /F cone along the direction of the 
photon beam, where F is the bulk Lorentz factor of 
the forward shock, but enhance the emission at larger 
angles, which leads to a delayed arriving time of the 
IC photons relative to the flare photons (Wang et al. 
2006; Fan et al. 2008). On the other hand, the self 
inverse-Compton emission of X-ray flares may give 
rise to a high-energy component which should have 
similar light curves as the flares (Wang et al. 2006; 
Fan et al. 2008; Yu & Dai 2009). During the early af- 
terglows, the afterglow synchrotron emission can also 
produce an extended high-energy component, which is 
thought to be responsible for the decaying long-lived 
high-energy emission detected from several GRBs by 
FermifLAI (e.g. Kumar & Barniol Duran 2009, 2010; 
GhiseUini et al. 2010; Wang et al. 2010; De Pasquale 
et al. 2010; Corsi et al. 2010; He et al. 2011, Liu & 
Wang 2011). 

GRB 100728 A is the second case (after GRB 
090510) with simultaneous detections by Swift and 
FermifLPC^ (Abdo et al. 2011). In this paper, we first 
analyze the Swift XRT data and the FermilhPC^ data of 
GRB 100728 A, and obtain the X-ray and high-energy 
(>100 MeV) Hght curves (§2). The early (167-854 
seconds) X-ray afterglow exhibits intense and long- 
lasting flaring activity, with a total of 8 flares (Abdo 
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et al. 2011). Two X-ray flares are also seen at the 
late stage (80-167 seconds) of the prompt phase (see 
Fig.l), indicating that the X-ray activity continues 
from the prompt phase to the afterglow phase. Inter- 
estingly fermZ/LAT detected high-energy gamma-ray 
emission during the X-ray flaring activities (Abdo et 
al. 2011). These simultaneous X-ray and GeV obser- 
vations offer a good case to study the nature of high 
energy emission. In §3, we confront the observations 
with the theoretical models and find that the external 
inverse-Compton (EIC) scattering of X-ray flare pho- 
tons by electrons in the forward shock provides the 
best explanation for the GeV emission. We calculate 
the light curve of the EIC emission by taking into ac- 
count the anisotropic scattering effect, and compare it 
with that of the observed GeV emission by FermifLAT 
in §4. Finally, we give our conclusions in §5. 

2. Observational facts of GRB100728A 

Bright GRB 100728 A was triggered by both 
Swift/BKr and Femii/GBM with Tgo 163 s. Intense 
X-ray flares are observed by Swift/XRT and significant 
GeV photons are detected by Fermi/LAT in the mean- 
while during the early afterglow phase. Our study 
mainly focuses on this phase to model both the X-ray 
and GeV behaviors. The XRT data are reprocessed 
by our own codes (Zhang, B.-B et al. 2007). Pass 7 
LAT data are retrieved from the Fermi LAT data server 
and reprocessed using a likelihood method. We se- 
lected "transient"(evclass=0) LAT photons during the 
prompt emission phase and "source" (evclass=2) LAT 
photons during the afterglow phase in a 20-degree 
circular region. The time bins are judged according 
to the separation among the prompt emission, flare 
and pure afterglow phases. Our results are generally 
consistent with Abdo et al. (2011). The observa- 
tional facts related to our modeling in this work are 
summarized as follows: i) The burst is bright with 
a fluence of ~ 1.3 x lO-^'ergcm"- in 10-1000 keV; 
ii) Ten successive flares in XRT (from ^ 80 s to 

854 s) are detected, the time-average spectrum of 
these flares can be described by a Band function with 
a = -1 .06 ± 0. 1 1, ^ = -2.24 ± 0.02 and a peak energy 
Spk = 1.0!|] jkeV (Abdo et al. 2011); iii) The spectrum 
of the GeV emission is hard with a photon index of 
Flat =-1-4 ±0.2 (Abdo et al. 2011); iv) The flux 
of GeV emission during the time t ^ 167-854 s is 
-Plat ^ (5.8±4.5) x lO'^^ergcm'^s"'. In what follows, 
we investigate which theoretical scenarios can explain 
these observations. 



3. Confronting the observations with various 
models 

3.1. Afterglow synchrotron emission scenario 

The underlying X-ray afterglow flux during the 
flare period should be Fx.af < 10~^ergcm"^s"' at time 
f ^ 500 s, according to Fig.l. The post-flare X-ray af- 
terglow has a decay index a2= 1 .07 ± 0.05 before t 
lOks and an average photon index of Fx = -2.07 ±0.09 
(i.e., (32 = -1.07 ±0.09 for the spectral index in the 
convention /v t^-v^-) (Abdo et al. 2011, Evans & 
Cannizzo 2010). Since the temporal index and spec- 
tral index satisfy a2 — ^^j^, i.e., the closure relation 
for the afterglow emission in the fast-cooling regim^H 
for constant-density interstellar medium (ISM) case 
(Zhang & Meszaros 2004). The power law index of 
electron number distribution p ( i.e., drie/d^e oc 7J''), 
can then be obtained p = -2/^2 — 2.2 (Sari et al. 1998, 
Zhang & Meszaros 2004). Thus, extrapolating it to 
the LAT energies, the GeV flux from the forward 
shock should be Fcev.af = Fx.aKlGeV/lkeV)-^''^^ < 
0.6 X 10~^ergcm"^s"'. This flux is about one order 
of magnitude lower than the observed GeV flux, disfa- 
voring the forward shock synchrotron emission model. 
Moreover, the hard spectrum of the GeV emission 
(Flat = -1 .4 ± 0.2, much harder than the X-ray spec- 
trum with Fx = -2.07 ± 0.09) cannot be explained by 
this model. 

3.2. X-ray flare self-IC scenario 

As argued in Abdo et al. (2011), the LAT emis- 
sion extends over the flaring period in the early after- 
glow phase, rather than mainly originated during the 
higher-significance flares. A cross-correlation analysis 
between the LAT (diffuse-class) and XRT light curves 
does not detect any significant temporal correlation or 
anti-correlation between the two data sets (Abdo et al. 
2011). These facts, if true, would disfavor the self-IC 
scenario for the GeV emission, since this scenario pre- 
dicts a tight temporal correlation between X-ray flares 
and GeV flares. However, we should note that the GeV 
emission signal is not sufficiently strong to allow one 
to draw a firm conclusion. 

The hard GeV spectral shape can be accounted for 
by a self-IC component only if it peaks in or above the 
LAT energy window, i.e. Spjc ^ IGeV. As the peak 
energy of the IC emission and the seed photon emis- 



Fast-cooling means that the electrons producing the observed radia- 
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sion are related by epjc ~ 27^ep.syn ^ '^ll.e^k with the 
observed peak energy Cpk = 1 .OkeV (Abdo et al. 201 1), 
this requires > 700, where 7m is the characteristic 
Lorentz factor of the electrons. Since the characteris- 
tic Lorentz factor of electrons (in the comoving frame) 
is 7m = 1.8 X 103(p-2)/(p- l)e,(r,h- 1) (Wang et al. 
2006), where is the equipartition factor of electrons 
and Fsh is the Lorentz factor of the internal shock, this 
would imply e, > O.Mp-l)/(p-2)(T,t- l)"^ Since 
the relative shock Lorentz factor may be of order unity 
(Wang et al. 2006), i.e. Fsh ^1-2, the inferred value 
of Ee is > 1 for /? ~ 2-3, which is too large to be re- 
alistic. We note that this scenario can not be excluded 
if the Lorentz factors (Fsh) of the internal shocks for 
these X-ray flares are very large. 

3.3. The external IC scenario 

The hard spectrum can be more easily accounted 
for in an external IC (EIC) scenario, in which the flare 
photons are scatterred by external forward shock elec- 
trons. This is because the Lorentz factors of the elec- 
trons in the forward shock are larger by a factor of 
F (i.e., the bulk Lorentz factor of the forward shock), 
compared to the internal shock case, and hence the EIC 
emission peaks at much higher energies. In the adia- 
batic case, the shock Lorentz factor can be derived by 



(1) 



where we use the convention Qx = 2/10^ in cgs units 
hereafter (Sari et al. 1998). Adopting equation ([T]), we 
can obtain the typical Lorentz factor of the post-shock 
electrons at time f as 7^ ~ 1.8 x 10-^{p-2)/{p-l)eeT = 

, E is 



1.3 X lO'^ee.-ifpEli^^nQ^^^tj^^^ , where fp 



6(p-2) 
ip-l) 



the kinetic energy of the blast wave and n is the num- 
ber density of the circum-burst medium. By inputting 
the value of the observed peak energy Epk = l.OkeV 
(Abdo et al. 2011) into hv^^ = 2-/l^epi, (Saii & Esin 
2001, Wang et al. 2006), we can obtain the peak of the 
observed EIC vF^ flux as 

/..^■^=4xl0^eV4,4/^«-'/V^' (2) 
for p = 2.2. 

The high flux of the X-ray flares will result in en- 
hanced cooling of the electrons in the forward shock. 
Adopting equation ([T]) and the shock radius as R = 
4F^cf(Waxman 1997), the energy density of X-ray 
flare photons in the forward-shock frame is 



U^ = dlF^/{T^R^c) 

nc!yii7-3/4 3/4„ ,2 A/4 _i _2 

= 0.34^53' Kg' Fx.-%dl2ih ^rgs 'cm^, 



(3) 



where Fx is the observed flare flux, di^ is the lumi- 
nosity distance of the burst. Considering the case in 
which EIC cooling is dominant, the cooling power 
is P(7,) = \aTcT^-ilUl^ (Sari et al. 1998), then the 
cooling Lorentz factor of electrons can be obtained by 
equation T^cmgC^ = P{ jc)t (Sari et al. 1998), i.e.. 



7c ~ 3mi,c^ / (16(TjcU^rt) 



: 5 X 10^£, 



2^.5/8^-5/8 „_i 



53 "0 



- X-8"L,28'3 



-7/8 



(4) 



Thus, adopting equation (|4]i and £pk = l.OkeV, the 
cooling break in the EIC spectrum can be obtained as 
(Sari & Esin 2001, Wang et al. 2006) 



hi^^'"- = 27c^epk 
= 5 X 10^eV£, 



5/4, 



,-5/4.-7/4 p_2 ,-4 
'3 ■'^X,-8"L,28- 



(5) 



Below i^rn the EIC spectrum has a photon index of 
Teic = -3/2 if i/f^'^ <v < I'm'^, or has the same index 
as the low-energy index (a = -1.06 ± 0.1 1) of the X- 
ray flare spectrum if < mm{v^^'~^ ,v^'~^) (Sari & Esin 
2001, Wang et al. 2006). Therefore, below y^'^, the 
photon index of the EIC emission is consistent with the 
hard spectrum of the observed GeV emission. Flat = 
-1.4±0.2. 

Using R = 4F^cf for the radius of the forward shock 
and the spectral peak epk = 1 .OkeV (Abdo et al. 201 1), 
we can derive the peak spectral flux of the EIC emis- 
sion as (Sari & Esin 2001) 

ff^ = r,..{^y = \o,nR[^)K 

= 2 X 10-^^fe„£53/'*no'''*f3''Vx-8ergcm-2s-' Hz"', 

(6) 

where Tssc = \<7T:nR is the scattering optical depth of 
the forward shock shell, aj is the Thompson cross sec- 
tion, ka is the correction factor accounting for the sup- 
pression of the IC flux due to the anisotropic scattering 
effect compared to the isotropic scattering case. It is 
found that this correction is mild with kg ~ 0.4 (Fan & 
Piran 2006; He et al. 2009). 

Therefore, the peak flux of the EIC emission for 
the fast cooling case (i.e., v^^'^ < i^^'^) is (Sari & Esin 
2001) 



„EIC rEIC 



ft. 



. ..EICrEIC { 
■ '^m Jp \^i/Eic J 



-1/2 



= 7 X 10 ^^^kaE^iee-itj^ d[\gergs 'cm 



(7) 



From the spectrum of the observed GeV emission 
at f = lO-'s, we have two constraints, i.e.. 



(8) 
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and 

''m'^fm^ ^ 5 X lO-'^ergcm-^s"' . (9) 

Inputting equations ^ and (|7]i into the above two con- 
straints, respectively, one can get the following con- 
straints 

nT (10) 

and 

0.01k-' E^Ul^^. (11) 

The above constraints are consistent with the post- 
flare X-ray afterglow observations. We attribute the 
post-flare X-ray afterglow to the synchrotron forward 
shock emission. Since its temporal index and spectral 
index satisfy the closure relation a2 = as men- 

tioned in section 3.1, the X-ray afterglow emission is 
in the fast cooling regime. According to equation (9) in 
Zhang et al. (2007), the synchrotron X-ray afterglow 
flux is 

= 2 X lO-^^ergs-i cm'^ 

(12) 

where gp = 8.6 (^)'' (33 x lO'^)"^, and the IC 

parameter is F = (-1 + ■\/l+4^Q/eg)/2 with rj being 
the radiation efficiency (Sari & Esin 2001). According 
to Sari et al. (1998), the afterglow cooling Lorentz fac- 
tor of electrons is 7f = 1.6 x \O^E';l'^nQ''^e-^'_/J^{\ + 
7)"'. If the IC parameter Y is not too large, we 
have > 7„, at f ~ lO^s, implying a slow cooling 
case, where the radiation efficiency is 77 = (7^^*/ 7,,,)^ '^ = 
OMpE^;''nlr''^\0,e':yt''^^\l+Yy-^ with = 
l.l^-^-P (Sari & Esin 2001). At t = 10"^ s, for Dl = 
2 X lO^^cm (i.e., at redshift z=l), with ~ 0.3, - 
0.1, n = Icm""* and £ ^ 3 x lO^-'erg, the IC parameter 
is F 2± 1. For the above parameters, the flux derived 
from equation (fT2] i is consistent with the observed X- 
ray flux vF^^^ ^ 3 x 10""ergcm"-s"'. 

Note that in the above estimate we have implicitly 
assumed that the X-ray flare flux is sufficiently strong 
that the electrons in the forward shocks are cooled 
down by the flare photons (i.e. > 7c). Neverthe- 
less, in the regime that electrons cool slowly ( i.e., 
7m < 7c in the slow cooling case), one can also have 
an EIC spectrum as hard as the low-energy spectrum 
of the X-ray flare (i.e., Teic ^ a = -1.06 ±0.11) if 
hi/^'^ > IGeV. Such a situation will be included in 
the numerical modelings in §4. 



3.4. The afterglow synchrotron self-Compton 
emission scenario 

Predictions were made that the synchrotron self- 
Compton (SSC) emission from the forward shock 
electrons can produce GeV emission during the early 
afterglow in some parameter space (e.g. Zhang & 
Meszaros 2001). In the case of GRB 100728A, how- 
ever, the illumination by X-ray flare photons enhances 
the cooling of the forward-shock electrons, which in 
turn suppresses the afterglow synchrotron and SSC 
emission. We expect that the ratio of the EIC lumi- 
nosity to the SSC luminosity is Leic/i'SSC = f^x/^^syn' 
where and U^^^ are, respectively, the energy den- 
sity of flare photons and synchrotron photons in the 
comoving frame of the forward shocks. During the 
flare activity period, because the flux of flare photons 
is significantly larger than that of the synchrotron X- 
ray afterglow, i.e. f/^ > ^^syn' "^SC emission should 
be subdominant compared to the EIC emission. 

4. Numerical modeUng of the GeV emission 

In this section we calculate numerically the flux 
of the EIC emission of the x-ray flare photons scat- 
tered by forward shock electrons, taking into account 
the anisotropic scattering effect. For a photon beam 
penetrating into the shock region where electrons are 
isotropically distributed, the EIC emissivity of radia- 
tion scattered at an angle 6'sc relative to the direction 
of the photon beam in the shock comoving frame is 
(Aharonian & Atoyanl981, Brunetti 2000, Fan et. al. 
2008, He et al. 2009): 

e--(.',cos^sc) « ^ i:;r d^^w. !<:: I? 

(13) 

where ly'^ and v' are respectively the comoving- 
frame frequencies of the photons before scatter- 
ing and after scattering, is the Lorentz factor of 
scattering electrons, ^ = hv' /(jenieC^), bg = 2(1 - 
cos9^c)lehv'^l (fneC^), and hv'^ <^ hv' <^ ^enieC^bg + 
bg). Integrating equation (flJl l over the angle 0sc for 
the whole solid angle, one can reduce equation (T% 
to the equation for the case of isotropically distributed 
electrons and photons. The maximum Lorentz factor 
of electrons 7inax is obtained by equating the electron 
acceleration timescale to the cooling (including syn- 
chrotron cooling and IC cooling) timescale. fj^ is 
the flux density of seed photons at the frequency 1/'^, 
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which are X-ray flare photons in our case. The lowest 
and highest frequencies of seed photons in the shock 
frame can be calculated via = 0.3keV/r and 
'^smax - lOkeV/r, according to the observation range 
of Swift XRT, i.e., 0.3keV- lOkeV. The light curves 
of these X-ray flares are modeled by power-law rises 
and decays in the calculation, as shown by the green 
Unes in Fig.l and Fig.2. The spectra of these X-ray 
flares are modeled by the Band function as given in 
section 2. 

The electron distribution dtie/d^te in Eq.(13) de- 
pends on where the cooling break 7c is. When the 
circum-burst density is sufficiently low, the forward 
shock electrons suffer from weak cooling by X-ray 
flares, so we have 7m < 7c. Otherwise, we have 
7m > 7c- The distribution of electrons is given by (Sari 
et al. 1998) 



for 7c < 7m < 7max, and 



oc 



le ' 



7c < 7e < 7m 
7m < 7e < 7max 



7m < 7e < 7c 
7c < 7e < 7max 



(14) 



(15) 



for 7m < 7c ^ 7max. 

The flux density of the EIC emission at a frequency 
u is given by (Huang et al. 2000, Yu et al. 2007) 



Jo 



9i ^'EIC 



s'^^^iv/D,cose)D^ 



InR^sme COS ed9, 



(16) 

where D = l/[r(l-/3cos6')] is the Doppler factor, (3 
is the velocity of the forward shock, 9 is the angle be- 
tween the motion of the emitting material and the jet 
axis, and dj is the opening angle of the jet. Since 9 also 
represents the angle between the injecting photons and 
the observed scattered photons in the observer frame, 
it is related with the angle 9sc in the comoving frame 
by cos ^sc = (cos6'-/3)/(l-/3cos0). 0jisa function of 
time with an initial value of 6jfl. Since the integration 
of Eq. (13) is performed over the equal arrival time 
surface (EATS), the radius R of the shock is a function 
of 6, which is determined by 



■ B cos 9 



13c 



-dR = constant 



(17) 



within the jet boundaries (e.g. Waxman 1997; Granot 
et al. 1999). 

The calculated light curves and spectra of the EIC 
emission are shown in Fig.l and Fig.2 for high (n = 



Icm"-') and low (n = O.Olcm"-') circumburst density 
cases respectively, in which the electron distribution is, 
correspondingly, in the fast-cooling and slow-cooling 
regime. In both cases, the EIC emission is the domi- 
nant component at GeV energies during the X-ray flare 
activity period. In the high density case (Fig.l), the 
electrons are in the fast cooling regime during the X- 



ray flare period and we have I'c < '^at < i^m > lead- 
ing to a hard GeV emission spectrum with a photon 
index Fcev = -1-5, as shown by the inset spectra of 
Fig. 1. In contrast to the high temporal variabiUty of 
the x-ray flares, the EIC emission is very smooth. It is 
clearly seen that the EIC emission starts later than the 
earhest flares, which can explain the non-detection of 
GeV emission during the earliest two flares (80-167s). 
The EIC emission continues after the X-ray flare ac- 
tivity shuts off, which may explain the marginal detec- 
tion of GeV emission after the X-ray flare activity, as 
reported in Abdo et al. (201 1). 

Fig.2 shows the case of a low circumburst density 
with n = O.Olcm"^. In this case, a large shock radius 
results in a lower flare photon density in the shock co- 
moving frame so that the EIC cooling is not impor- 
tant. In this case, we will have 7m < 7e, and thus 
ly^Aj < ^m *^ < t'^'^- The spectrum of the EIC emission 
below v^'" has the same slope as that of the seed pho- 
tons, i.e. the photon index is similar to the low-energy 
photon index a of the Band function describing the 
spectrum of X-ray flares. Such a spectrum is also con- 
sistent with the measured spectrum of the GeV emis- 
sion within the error bars, as shown in the inset plot 
of Fig.2. In the low density case, the afterglow syn- 
chrotron emission dominates over the SSC emission in 
the GeV band in the early time. 

5. Conclusions 

For the first time, GeV emission from a GRB dur- 
ing the X-ray flaring activity has been detected by the 
Fermi/L AT, in GRB 100728 A. The temporal coinci- 
dence between the GeV emission and the X-ray flares 
suggests that the GeV emission should be related to the 
flares in some way. Here we have shown that an EIC 
scenario, where X-ray flare photons are up-scattered 
by electrons in the external forward shocks, provides 
the best explanation for the GeV emission, supporting 
the earlier prediction of GeV emission from GRB X- 
ray flares (Wang et al. 2006). The hard spectrum of 
the GeV emission can be readily explained by the EIC 
emission. The delayed behavior of the GeV emission 
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relative to the X-ray flares also naturally arises in the 
EIC scenario. Synergistic observations between Swift 
and Fermi should be able in the future to find other 
similar cases, allowing further tests of this EIC sce- 
nario. 
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Fig. 1. — Light curves of the GeV emission and X-ray afterglow of GRB 100728 A. The filled black squares represent 
the observed X-ray flux by Swift XRT (0.3-lOkeV) and the empty black squares represent the high energy data and 
upper limits from Fermi/LAT (100MeV-30GeV). The green solid line represents the input light curve of X-ray flares 
that we used in the calculation. The blue solid line represents the synchrotron afterglow emission in 0.3-10 keV, while 
the red solid line represents the high energy emission (in lOOMeV- 30Ge V) produced by the external inverse-Compton 
scattering of the X-ray flare photons. The blue dotted line and purple dash-dotted line represent, respectively, the 
afterglow synchrotron emission and the afterglow SSC emission in 100MeV-30GeV. The black solid line is the sum 
of the three high energy components (represented by the red solid line, blue dotted line and purple dash-dotted line, 
respectively). For clearness, we re-scaled the flux of all the high-energy components by multiplying a factor of 0.01. 
The inset plot shows the spectra of the three high-energy components and the LAT data taken from Abdo et al. (201 1) 
during the X-ray flaring period (167-854 s). £ = 3 x 10^^erg,n = lcm-\p = 2.2, T = 250,6*^,0 = 0.1, e, = 0.1, = 0.04 
and z = 1 are used in the calculation. 
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Fig. 2. — The same as Figure 1, but for a low-density case (n = O.Olcm"-'). £ = 1 x 10^^erg,n = O.Olcm"^,/? = 2.2,r = 
400,9 jfl = 0.07, ee = 0.2, eB = 0.1 and z = 1 are used in the calculation. 



7 



